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Abstract 

The global stability of i?-charged AdS black holes in a grand canonical ensemble is 
examined by eliminating the constraints from the action, but without solving the equa- 
tions of motion, thereby constructing the reduced action of the system. The metasta- 
bility of the system is found to set in at a critical value of the chemical potential which 
is conjugate to the i?-charge. The relation among the small black hole, large black hole 
and the instability is discussed. The result is consistent with the metastability found in 
the AdS/CFT-conjectured dual field theory. The "renormalized" temperature of AdS 
black holes, which has been rather ad hoc, is suggested to be the boundary temperature 
in the sense of AdS/CFT correspondence. As a byproduct of the analysis, we find a 
more general solution of the theory and its properties are briefly discussed. 
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1 Introduction 



The global structure of a black hole action is the interest of this paper. Examination of 
this, in general, is a formidable task and the difficulties include the large number of the 
fields involved and the complicated nature of the action itself. For instance, in a single 
scalar theory at classical level, one simply plots the potential of the theory with respect to 
the expectation value of the scalar field and examines the global minimum of the theory. 
Similar program for a theory that involves gravity is not as straightforward due to the 
complications. Therefore, a black hole action is commonly studied by evaluating the values 
of it on known black hole solutions, by this means constructing the on-shell action of the 
theory. The on-shell action is valid in an infinitesimal neighborhood of the black hole saddle 
point and the local behavior of black holes, such as the first law of thermodynamics and 
local thermodynamic stability, can be studied from the action. This method, however, is 
not adequate for examining the global structure of a black hole action. 

In Reference pQ, Braden et al. devised the technique that yields a simplified off-shell 
action, called "reduced action" of a black hole system. The procedure for constructing the 
reduced action is as follows. One starts with the action of a system and place it in a finite 
box. Then, restrict the analysis to a class of geometries that satisfy certain conditions, such 
as static spherical symmetry. The restriction results in a specific form of the metric and this 
is analogous to a metric ansatz that one makes in solving the Einstein equation. In addition 
to the geometrical restriction, one imposes the black hole "cigar" topology on the spacetime 
manifold, giving a condition that the metric must satisfy. Next, the thermodynamic data 
is encoded on the wall of the box. The data, in the case of a grand canonical ensemble, is 
the temperature and chemical potential. (We will show exactly how to encode the data in 
Section f3 . 2 [) . As will be explained in Section [3. 3\ the physical states of the system in general 
obeys constraint equations including the Hamiltonian constraint, and they are solved and 
eliminated from the action. Finally, the action functional is integrated with the topology 
condition and the boundary data, resulting in the simplified reduced action. Notice that 
the equations of motion have not been used in reducing the action. Thus the reduced action 
is valid throughout the physical parameter space and a black hole solution corresponds to 
a saddle point of the action in the parameter space. Therefore, this action captures the 
global structure of the theory. 

After reviewing the five dimensional i?-charged black hole solution and its properties in 
Section [21 in Section El we apply the procedure sketched in the previous paragraph to this 
black hole system in a grand canonical ensemble. This black hole system can be obtained 
through 5 5 compactification of ten dimensional Type IIB supergravity with gauged 50(6) 
global symmetry [21 [3] . A special case with decoupled scalars can be obtained rather simply 
by "twisting" the five sphere upon the compactification, as demonstrated by Chamblin 
et al. in Reference It has also been shown by Cvetic and Gubser in Reference [5] 
that the black hole geometry arises from the near horizon limit of rotating D3-branes in 
ten dimensions with the rotation in their transverse directions. Therefore this system has 
drawn numerous interests in connection to string theory and the examination of the global 
structure is worthwhile in its own right. However, there is an additional and more specific 
motivation for studying the global stability of the theory. 

According to the AdS/CFT correspondence [El El [8], the i?-charged black hole system in 
a grand canonical ensemble is conjectured to be the dual of strongly coupled thermal N = 4 
supersymmetric Yang-Mills theory in the large -/V c limit with the chemical potentials corre- 
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sponding to the maximal Abelian subgroup of the SU(4) i?-summery. The weak coupling 
analysis of the field theory in Reference [9] shows that, despite the difference in the coupling 
regime, the phase diagram in the parameter space of temperature and chemical potential 
closely resembles the one obtained in the gravity analysis of Reference [5]. However, in 
the phase diagram of the field theory, the additional structure has been found where the 
equilibrium state of the system becomes metastable beyond a certain critical value of the 
chemical potential. 

The metastability in the field theory was found by examining the global behavior of the 
scalar field effective potential at one-loop level. When the system is in the deconfined plasma 
phase with a sufficiently low value of the chemical potential, the origin of the scalar field 
space is the absolute minimum of the (one-loop) effective potential and defines the stable 
ground state of the theory. When the chemical potential is raised beyond the critical value, 
the effective potential goes down to minus infinity at the asymptotic region of the scalar 
expectation value. The minimum at the origin persists to exist in this phase but it now 
only is a local minimum with a potential barrier separating it from the unstable direction. 
The state at the origin, therefore, is metastable in this phase and the true ground state of 
the theory does not exist. Further rise in the chemical potential would eventually lead to 
the local instability of the state where the barrier that separates the local minimum and 
the unstable direction disappears. 

We suspect that the metastability is caused by the violation of the BPS inequality 
expected for the states of the field theory. The chiral primary fields of M = 4 super- 
Yang-Mills theory contain the half-BPS operators that are characterized by their scaling 
dimension equaling their 72-charges and in general, the dimension of an operator must be 
larger than or equal to the charge. When the theory is conformally mapped to a sphere, 
the operator and the dimension are mapped to a state of the theory and the energy of 
the state, respectively. Therefore, the half-BPS states have the energy being equal to the 
i?-charges. When we introduce the chemical potential conjugate to the charge, we expect 
that the product of the charge and the chemical potential contributes to the total energy 
of the state. (Such an argument has been made in Reference [10] and further generalized in 
Hence if the chemical potential is raised higher than one (in the units of the radius of 
the sphere), we expect that the BPS inequality is violated and something should go wrong 
with the theory. The metastability of the field theory analysis sets in exactly at this critical 
value of the chemical potential. 

Since the dimension of the BPS operator is protected against the quantum corrections 
and change in the coupling constant, we expect similar behavior in the gravitational analysis 
which is supposed to be the strong coupling region of the same theory. But it is clear that 
the local on-shell analysis that has been carried out in Reference [5j is not able to detect 
the metastability of the black hole system. This is why we embark on the program outlined 
above. Analyzing the resulting reduced action, we will find the metastability in Section SJ 

In Section [U a few issues that come up during the analysis are discussed. First in 
Section \5. 11 the mysterious "renormalized" temperature that has been adopted in literature 
on the AdS black hole systems is suggested to be naturally understood in the context of the 
AdS/CFT correspondence. Then Section \5 . 2 1 presents a more general solution of the theory 
than what has been known. This solution is obtained as a byproduct of the main analysis 
and its properties are briefly discussed. Section 15.31 addresses a subtle issue of the physical 
interpretation of the counterterm that does not diverge as the wall of the box is pushed to 
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infinity. Finally in Section 15.41 we contrast the revised phase diagram of the gravity side to 
the one in the field theory side. 

2 Review of i?-Charged Black Holes 

We start by reviewing the i2-charged black hole solution and their thermodynamics, mainly 
following Behrndt et al. [12j and Cvetic and Gubser [5]. The analysis on the instability of 
the black holes is slightly extended from the work of Cvetic and Gubser [5] for the cases 
with different charge configurations 

The -R-charged black hole of our interest is the solution to five dimensional M = 8 
gauged supergravity where the gauged symmetry is the U(l) 3 truncation of the 50(6) 
global symmetry that arises from the 5 5 compactification of ten dimensional supergravity. 
The solution is given by the metric 

ds 2 = -H(r)- 2 / 3 f(r)dt 2 + tf(r) 1 / 3 [/(r^dr 2 + r 2 dn 2 ] , (1) 

where d^l 2 is a metric on the unit three sphere and the function H (r) is 

H(r) = H 1 (r)H 2 (r)H 3 (r) with Hj{r) = 1 + ^ , (2) 

and 

fir) = 1 - f 2 + ^H(r) . (3) 

The physical interpretations of the parameters qj, ro and R are as follows. As mentioned 
in the introduction, it is possible to interpret this solution as the near horizon limit of 
rotating D3-branes in ten dimensions and in this perspective, the parameters qj with the 
dimension length squared are related to the three possible transverse angular momenta of 
the D3-branes. (See Reference |5j.) The parameter ro in the function / measures the "non- 
extremality" of the black hole in that ro = corresponds to the extremal BPS-saturated 
solution. The non-extremality parameter can be expressed in terms of the horizon radius 
th by solving the equation f(rn) = 0, explicitly, 

rl = r 2 H {l+ r f 2 H{r H) y (4) 

The parameter R is the length scale of the theory. Note that the solution approaches 
to the AdS metric at the r — > oo asymptotic region and we see that the parameter R is 
the curvature radius of this asymptotic AdS space. Therefore, it is also related to the 
cosmological constant. 

This theory contains scalar fields, Xj(r) with j = 1, • • • ,3, satisfying the constrain^ 

Xi(r)X 2 (r)X 3 (r) = 1 . (5) 

The extended results were also reported in [9]. 
2 As in Reference [12] the location of the index j, either up or down, does matter sometimes. In this 
work we exclusively use the lower index, which corresponds to the upper index of the reference, purely for 
the presentation reason. The notation is consistent within this work, however, if the reader is interested in 
comparing expressions with the reference, the convention difference must be watched. 
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The solutions to the scalar fields are 

X 3 (r)=H 3 (r)- 1 H{r)^ . (6) 

This theory contains three dimensionless Abelian gauge fields arising from the gauged 
i7(l) 3 and the solution takes the form 

A j Ar) = Sv,a(^i 2^—) > W 

\ r H + 1j r 2 + qj J 

where v is a spacetime index, the parameters ej with the dimension length squared are the 
electric charges of the system under the Abelian fields, the index j runs, again, from 1 to 
3 and we have set Aj u (rn) = for the reference point of the potentials. Observe that the 
time-component of the gauge fields at the coordinate r is the electric voltage between the 
horizon and the location r. We denote the potentials at the boundary with r = oo by \ij 
and they are 

H = 4 > (8) 

Rr 2 H + Qj 

where the length scale R is inserted for the dimensional reason. The parameters qj and ej 
are related through 

r 2 

e] = qj (r 2 H + qj ) l + ^§ Y[H k (r H ) 
1 Mi 

where the repeated index j is not summed over. 
The ADM mass M of the black hole is 



(9) 



v 7 = 1 ' 



(10) 



r- 

.2 



where LU3 = 2tt is the volume of the unit three sphere and G5 is the five dimensional 
Newton's constant which has the dimension of length cubed. This ADM mass assumes 
the AdS background without the presence of the black hole, i.e., the mass M is the mass 
difference between the AdS black hole spacetime and the thermal AdS space without the 
black holeU The entropy S is a quarter of the horizon area as usual: 

Q _ A _2vrH7^)_ (n) 



4G 5 4G 5 
The inverse Hawking temperature (3jj is given by the expression 



2nR^ H ^H(^Hj 

In Reference [5], the thermodynamic properties of the black hole in the grand canonical 
ensemble is examined by introducing the Gibbs free energy 

/Gibbs =(3 h (m-J2 ejVj) - S , (13) 
^ 7=1 ' 



3 In Reference [13] . Buchel and Pando Zayas point out that the ADM mass actually is divergent due to 
the presence of the gauge fields. However, the Brown- York '14 a type computation by Liu and Sabra [15] who 
include the "finite counterterm" yields the mass as presented in Equation IjlOp . The finite counterterm is 
discussed further in Sectiorl5T3l 
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where the length scale R is set to 1. The parameters (3h and \ij are the fixed data of the 
ensemble and the black hole corresponds to the saddle point of this function. Notice that 
the potentials fj,j are taken to be the chemical potentials of the grand canonical ensemble. 

As mentioned earlier, the ADM mass is the difference between the black hole spacetime 
and the thermal AdS space without black holes. Thus, a change in the sign of the Gibbs free 
energy indicates that the ground state (or thermal equilibrium state) of the system switches 
from the black hole to the thermal AdS or vice versa. This is the well-known Hawking-Page 
phase transition [16J. The phase transition line of the phase diagram in (fJ>,T) parameter 
space can therefore be determined by the condition /Gibbs = 0. 

Cvetic and Gubser [5] further examined the local thermodynamic stability of the black 
hole solutions. At an instability point, the black hole solution that extremizes the Gibbs 
free energy becomes an unstable saddle point, as opposed to a local minimum. Therefore, 
when the Hessian of Icibbs with respect to r# and qj is evaluated at the black hole solution, 
the determinant of the Hessian vanishes at the critical local instability point. This condition 
establishes the instability line in the phase diagram. Even though this criterion does not 
determine which side of the instability line corresponds to the stable black hole, one can 
figure that by knowing that the region of high temperature and zero chemical potential 
must correspond to a stable black hole solution. 
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Figure 1: The phase diagram on the left is plotted for (qi,q2,qs) = (q, 0,0) and on the right is for 
(q,q,q). Both \i and T are measured in the units of 1/R, the inverse of the AdS curvature radius. The 
solid and dotted lines are the black hole instability lines whereas the dashed lines are the Hawking-Page 
phase transition lines. The regions between the dashed and dotted lines indicate the spinodal phases where 
the AdS space is energetically preferred but the black hole still remains locally stable. The temperature 
To = 1/tt on the left-hand diagram is where th = 0. 



The resulting phase diagrams are plotted in Figure [TJ We choose to plot the phase 
diagram in fi-T plane as oppose to Reference [5] whose choice is q-rn plane. We have 
also extended the analysis of the same reference to the cases with more than one q is 
turned on. The left-hand diagram is plotted for (qi,q2,Q3) = {l,0,0) and the right is for 
(?) The solid and dotted lines are the black hole instability lines whereas the dashed 
lines are the Hawking-Page phase transition lines. The regions between the dashed and 
dotted lines indicate the spinodal phases where the AdS space is energetically preferred but 
the black hole still remains locally stable. The general structure of the phase diagrams are 
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similar for the different charge configurations [including the case with (q, q, 0)]. However the 
low temperature regions are quite different. The temperature To in the left-hand diagram 
denotes the Hawking temperature at which rjj = 0. When we have three equal non-zero 
charges, To = 0, but for the case with the only one non-zero charge, To 7^ 0. The exact 
values of the Tq are listed in Table [1] for the different cases. Since the action is the difference 



(91,92,93) 


(9,0,0) 


(9,9,0) 


(9,9,9) 




1/vr 


1/27T 





Asymp. Slope 


tt/x/2 


7T 


2tt 



Table 1: The list of the values of To and the slops of the high temperature linear behavior for different 
cases. 

between the black hole and AdS solutions, the region of the phase diagram below To is 
completely undetermined. 

In all the cases, the instability lines and the Hawking-Phage phase transition lines merge 
at T = Tq and \i = 1/R where R is the AdS curvature radius. The solid instability lines 
rise linearly at temperature large compared to 1/R and the slopes of the linear behavior is 
also listed in Table [1] for the different cases. 

We mention that the instabilities of the black holes with q\ = q^ and qi = q2 = Q3 are 
obtained by assuming that the equalities are subject to the thermal fluctuations around the 
black hole saddle points. Technically, this means that the Hessian of Icibbs with respect to 
qj is computed by assuming they are the independent parameters, and then the Hessian is 
evaluated at the black hole saddle point with the relations q\ = q2 or </i = <?2 = 93- One 
may choose to consider the black hole with the constraint q\ = qi or q\ = q<i = q-$ which 
is not subject to the fluctuations. In this case, the action Icibbs is obtained by setting 
9i = 92 or q\ = qi = qz first, and then the Hessian is computed with respect to the single 
charge (and rn). The resulting thermodynamic system for this constrained case is stable 
everywhere and the instability lines are not observed. The constraints reduce the dimension 
of the parameter space and eliminate the direction of the instabilities. This is the reason 
why the instability found in Reference [5] [with the configuration (q, 0, 0)] is not observed 
in the analysis of Reference [4] [with the configuration (q,q,q)]. 

In concluding the review, we point out that the Gibbs free energy Icibbs is constructed 
out of the quantities of the black hole. The black hole is the solution of the equations of 
motion and thus Icibbs is made out of the on-shell quantities. Therefore, even though we 
expect the Gibbs free energy to be the action of the theory in an infinitesimal neighborhood 
of the black hole extremum, we have no guarantee that Icibbs has any meaning as the action 
of the theory away from the extremum. The function Icibbs is a local on-shell action in this 
sense and it is not useful in examining the global stability of the theory. The examination of 
the global stability requires the construction of the off-shell action that is valid away from 
the extremum. 

3 Construction 

3.1 The Action, the Metric, and the Topology 

We now embark on constructing the reduced action of the system. In this subsection, we 
first display the action and the metric ansatz to start. Then we assume the Euclidean black 
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hole topology on spacetime and deduce the resulting conditions on the metric. 



As shown in Reference [12] , the black hole solution reviewed in the previous section 
can be obtained from the following special case of the five dimensional J\f = 2 gauged 
super gravity action 



r 3 3 

/ = l — [ d 5 x^ -K - - V -^(d u X j )(di"X j ) - - V —-^Fa uv Ff v + 

8ttG 5 J v 2 2^2X? M 3 3 A^2X 2 JM 3 R 2 

.7=1 ■> .7 = 1 J 



(14) 



where only the bosonic sector is shown. We have adopted the Euclidean signature and 
the sign of the action is chosen so that the weight factor in the path integral appears as 
exp(— /). The Chern-Simons term is dropped from the action because we will be considering 
the electrostatic fields which have the form Aj v = Aj T 5 u fi, and the term identically vanishes 
for such fields. The determinant of the metric is denoted by g, and 1Z is the scalar curvature. 
As in the previous section, the fields Xj(x) are the scalars and they satisfy the constraint 

X 1 (x)X 2 (x)X 3 (x) = 1 . (15) 

The three gauged f(l) 3 C SO(6) Abelian fields Aj u have the corresponding field strengths 

Fj [iv • — Aj Aj , (16) 

where the indices /i and v are the spacetime indices and the semicolon denotes the geomet- 
rical covariant derivative with respect to a given metric, which, of course, can be replaced 
by normal derivatives in this case. The potential V is defined as 

v-=*t^- (») 

3=1 3 

Finally in the action, the parameter R is the length scale of the theory which is related to 
the cosmological constant of the asymptotic geometry. 

We consider the static gravitational system with spherical symmetry and thus take the 
metric ansatz 

ds 2 = b{rfdT 2 + tf(r) 1 / 3 [f{r)- l dr 2 + r 2 dn 2 3 ] , (18) 

where r is the periodic Euclidean time whose period is set to 2irR. Consequently, b(r) 
is a dimensionless function of r. The form of the metric is completely general under the 
static spherical symmetry but we choose the particular form that resembles the solution of 
the previous section in Equation ([1]). This choice is motivated by the fact that we want 
the resulting reduced action to reproduce the same form of the solution as in the previous 
section at the extremum. Without loss of generality, we may assume that the dimensionless 
function H(r) is positive definite and finite for finite values of r so that the horizon in the 
Minkowski signature is determined solely by the dimensionless function f(r). Notice that 
the static spherical symmetry requires the scalar fields Xj(x) be functions only of r as well. 

In the computations that follow, we place the thermodynamic system into a large box 
whose wall is located at r = re (the word "large" is used with respect to the length scale 
R). In the AdS/CFT language, this is the "cut-off boundary" of the system and will be 
removed by taking the limit rg — > oo. The lower end of the range of the parameter r is given 
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at r = rjj which is the horizon radius of the black hole in the Minkowski signature. We will 
loosely use the terminology, "horizon radius" , even if we are discussing the Euclidean black 
holes. 

In addition to the static spherical symmetry, we restrict our considerations to the action 
that describes the spacetimes with Euclidean black hole topology. The Euclidean black hole 
topology of the spacetime manifold is M 2 x S 3 and the boundary of the manifold has the 
topology S 1 x S 3 where the S 1 factor corresponds to the periodic Euclidean time circle. In 
the actual computation, we assume the topology S 1 x at the wall of the box, r = re, 
with the limit re — * oo in our mind. 

We first want to impose the topology M 2 x S 3 on the entire spacetime. However, the Euler 
characteristic of any odd dimensional manifold identically vanishes and it is not desired here. 
Thus, we fix one of the angular coordinates in and require this slice to have the topology 
M? x S 2 with the Euler characteristic \ equaling 2. This yields the condition^ 

fH-^KrH- 1 ^)'} 2 ] =0, (19) 

J r=ru 

where the prime denotes the derivative with respect to the radial coordinate r and the 
left-hand side of the equation is evaluated at r = th- We will see that this condition boils 
down to the expected condition for a black hole horizon, f(rjj) = 0. 

Next, consider the R 2 part of the spacetime. This is commonly called "the cigar" where 
the tip of the cigar corresponds to the black hole horizon and thus it is parametrized by 
the lower end of the radial coordinate r = rjj, whereas the other end of the cigar is the 
boundary S . In this topology, the S 1 time circle degenerates at the tip r = r# and thus 
we have 

b(r H ) = . (20) 

Also this topology requires the r-r slice of the spacetime to have the Euler characteristic 
X = 1 as the cigar is topologically a two dimensional ball with boundary. This condition 
parallels with the usually known "smoothness condition" at the tip of the cigar H Being 
cautious for the orientation of the S 1 time circle at the wall of the box, the Gauss-Bonnet 
theorem yields 

6 /jl/2^-l/6l = y R ( 21 ) 

J r=Tu 

There are terms evaluated at r = r#, one from the Gauss curvature part and the other 
from the geodesic curvature of the boundary curve, but they cancel and absent in the above 
expression. 



3.2 The Boundary Data 

We encode the boundary data of the system in this subsection. Those data will be the 
fixed parameters of the "histories" of the paths in the path integral and they will not be 
varied in the extremization of the action. Since we want to consider and fit the black 

4 The topological condition imposed here is the same as the one in Reference [T| and we have borrowed the 
result of the reference with a straightforward modification for the form of the metric that we have adopted. 

J Most commonly, the smoothness is achieved by setting the period of the S 1 time circle to a specific value. 
Though equivalent, we are fixing the period to 2irR and require the metric to satisfy a certain condition. 
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hole system in the context of the AdS/CFT correspondence, we have the metric satisfy the 
asymptotic AdS geometry near the boundary. Also we are interested in the grand canonical 
ensemble of the black hole thermodynamics, so the temperature and chemical potentials 
are the fixed parameters. If there is a physically equivalent field theory on the boundary, 
then those fixed black hole thermodynamic data should be identical to the ones in the field 
theory. Therefore, it is natural to encode the fixed parameters on the boundary where the 
geometry is fixed and the corresponding field theory is supposed to live. 

First, we want the asymptotic AdS geometry near the boundary. We impose the asymp- 
totic conditions on the metric ansatz (|18[) 



The parameter t is an arbitrary dimensionless positive constant. The motivation for intro- 
ducing the parameter I shall be given shortly in the discussion of the temperature below. 
One sees that with the analytic continuation, t = —i£r, the metric ansatz (|18p approaches 
to the AdS metric with the curvature radius R, up to the first order in tb- 

Next we discuss the temperature of the system. It is sometimes seen that an observer 
is not specified for the discussion of the black hole temperature. However, the inverse 
temperature is the circumference of the S 1 Euclidean time circle and it depends on the 
location of the observer. For our metric, the local, physical inverse temperature at r = r B 
is the proper length of the S 1 circumference 



Due to the asymptotic AdS condition on b(r) in (|22p . the physical temperature is red- 
shifted to zero as the wall is pushed to the boundary, i.e., as r# — > oo. This fact has been 
known since as early as Hawking and Page [16]. We could have discussed the thermody- 
namics defined in a finite box with r B < oo, but since we are interested in the AdS/CFT 
correspondence, we must take the limit r B — ► oo. 

In discussing the thermodynamics of AdS black holes in the infinite box, the unphysical 
temperature, so-called "renormalized temperature", has been used. (Among the works 
that adopt the "renormalized" temperature, we refer to one paper [17] in which this issue 
is clearly addressed.) The "renormalized" temperature is defined by the product of the 
physical temperature at r = r B and the red-shift factor of the black hole geometry so that 
the limit r B — ► oo yields a finite value. (Consequently, if we set the boundary temperature to 
the finite value, the temperature diverges everywhere in the bulk). Even though the physical 
meaning of the redefined temperature has remained unclear, there have been indications 
that this, at least formally, is the correct temperature to be used on the boundary. The 
indications include the correct forms of the heat capacity, first law of thermodynamics and 
the quarter area law of the entropy, all with the redefined temperature |16[ \T7\ 118] ■ In 
addition, our analysis below also requires a similar kind of redefinition in the temperature. 
Therefore, we introduce the new temperature (5 defined at the boundary by 



b{r B ) £ 



-5- , f{r B ) -> -Ms , H(r B ) 1 , as r B -> 00 . 



(22) 




(23) 




(24) 
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Here we see that the temperature is directly proportional to the arbitrary parameter I. 
Thus the temperature is a thermodynamic quantity that we can tune, not the system, as a 
grand canonical ensemble should be. We will discuss more on the redefined temperature in 
Section 15,11 



Finally, we encode the chemical potentials that are conjugate to the -R-charges of the 
system. In a grand partition function, a chemical potential couples to the conjugate charge 
which is the time component of a conserved current. Therefore, the chemical potential is 
directly related to the time component of the gauge field that minimally couples to the 
conserved current. Thus we define the chemical potentials at the wall r = r# as 



pij : = ilnexp 



; j> Aj /r 



(25) 



where the integral on the right-hand side is defined on the S 1 time circle at r = re and Aj 
are the Hermitian gauge field one-forms. In the definition, the length scale R is inserted for 
the dimensional reason. We choose the (Wick rotated) gauge field of the form 

A jv = A jT (r)6 v , . (26) 

This form assumes the electrostatic spherical symmetry. Then we have 

j> Aj = j>{-ilAj t )(idt/l) = <j> A jT dT = 2TrRA jr (r B ) , (27) 

where the analytic continuation t = —ilr has been done and Aj r now is pure imaginary. 
Therefore, we get 

2vri . . 
fij = —~—Aj T (r£i) . (28) 

P 

Now if we assume that the gauge fields Aj V {rB) are finite as re — * ooH then the 
chemical potentials fij at the boundary vanish because $ diverges due to the asymptotic 
AdS condition (|22p . Thus, similar to the case of the temperature, we define the boundary 
chemical potentials by 

^-rfe^K (29) 

This definition will be discussed further in Section 15.11 along with the temperature (3 defined 
in Equation (|2"3]). 

We add that since we have the condition (I20j) . we need to have the conditions 



A jT (r H ) = 0, (30) 

so that the proper local value of the gauge potential Aj t {vh) /b{rn) can be set to zero as 
the reference point of the potential. 



6 Alternatively, we may assume Aj ^(rs) are singular so that jlj remain finite. This can also lead to the 
same conclusion at the end. 
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3.3 The Constraints 

Our system has two sets of constraints; Gauss' law and the Hamiltonian constraint. The 
physical states of the theory must satisfy the constraints and we treat them fixed and 
not subject to the variation. The way they arise is very similar^ Among the Maxwell and 
Einstein equations, the constraints correspond to the non-dynamical equations, i.e., the ones 
without second order time derivatives. The constraints can also be viewed as the resulting 
equations of the variation with respect to Lagrange multipliers. In the electrodynamics, 
the Lagrange multiplier is the time component of the gauge field. In the ADM formalism 
|221I19| of the geometro dynamics, the shift and lapse functions are the Lagrange multipliers. 
Below, we bring the action (|14p to the form where the Lagrange multipliers explicitly appear 
and eliminate the constraints from the action, without imposing the dynamical equations 
of motion. 



We first consider the gravitational part of the action (|14D , 



I °—*kl*'*\ 11 - (31) 



One can verify that with the metric of the form (I18h . the scalar curvature 1Z can be written 

as 

4 K =6^73( 6V/1/2jfl/3 )' + GT - (32) 

where the last term is the ( T T )-component of the Einstein tensor and it explicitly has the 
form 



^^(fH'^iirH 1 ^)'} 2 -lYI . (33) 



2r 3 H 1 / 2 (rH 1 / e )' L 

Notice that this expression is independent of the function b. Now Iq can be expressed as 



l G = _||_27r J R| J dr(7/r 3 / 1/2 # 1/3 )' + J drbr 3 r 1/2 H 2 / 3 G T T \ . 



(34) 



To proceed further, we need to introduce the Gibbons-Hawking term [23J to the gravi- 
tational part of the action Iq. In general, the Gibbons-Hawking term is required so that a 
gravitational action has the consistent variational principle with respect to general variation 
of the metric whose induced metric on the wall of the box is held fixedcl The Gibbons- 
Hawking term is not unique and one is allowed to include arbitrary term that depends only 
on the induced metric on the wall which can be immersed in a different spacetime. This 



7 See Wald's textbook |19] for a pedagogical exposition of the constraints. It is also instructive to 
observe how the Gauss law and the Hamiltonian constraint are identically treated in References [20] and 
[21] . respectively. 

One might be slightly puzzled by the relation between the dynamical Einstein equations and the Hamil- 
tonian constraint in a static spacetime manifold, because the constraint is the initial value condition for the 
Cauchy development problem. The Hamiltonian constraint in a static spacetime determines the geometry of 
space-like leaves of the spacetime foliation but it does not determine the way those leaves should be "stack 
up" in the time-like direction. In other words, the Hamiltonian constraint cannot determine the function 
b(r) in the metric (|18|l and the determination of the function requires the full Einstein equations. 

8 If one considers more specific kind of the metric variation, namely, the metric whose induced metric on 
the wall of the box is held fixed and the first derivatives of the metric on the wall are also fixed, then the 
Gibbons-Hawking surface term is absent. 
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arbitrary term is commonly used to set the reference background geometry of the action. 
However, the introduction of the background geometry is arbitrary and sometimes it is not 
possible to set a seemingly natural background. Thus, it is more desirable to construct our 
action independent of the background, therefore, we do not include the ambiguous term 
here and instead, adopt the AdS/CFT motivated background independent counterterms as 
discussed in the next subsection. The Gibbons-Hawking term without the subtraction term 
has the form 

Igh = [ d^x^K , (35) 

o-KLx 5 Jq M 

where dM is the wall at r = r B with the metric 



ds\ = b(r B ) 2 dT 2 + H{r B ) l/Z r 2 B dnl , 



(36) 



and 7 is the determinant with respect to this metric. The factor K is the trace of the 
extrinsic curvature 



K, 



1 ( \ 

2 l^A 1 ;^ -i 



where rin is a normal unit vector of the wall, dM, and is given by 



n,, 



0) iji/6//i/2 ) 0j o, ) . 



A straightforward computation yields 

u 3 



Igh 



8ttG 5 

We combine Equations (|34|) and 

U) 3 



2vri? 



r=r B 



Ig + Igh 



87rGs 



13 



. After some algebra, we obtain 
f l/2 H -l/e {r 3 H l/2y 



(37) 



(38) 



(39) 



r=r B 



+ 2irr 3 H y/H(r H ) 



87tG 5 J 



(40) 



where the topological condition (|21|) and the definition of (5 in (|23p have been used. 
We proceed to the Maxwell part of the action 



lEM = 8^gJ '^lE^^^ 



fJbV 



(41) 



We define the Maxwell field tensors as in Equation (|16p and with our metric (|18p and the 
gauge fields ([26]) . the only nonzero components are 



(42) 



where the components with fj,,v = 1 are the radial r-component. We introduce auxiliary 
fields Ej(r) such that 



Iem 



UJ3 

8vrG, 



-,R [ drb-^f^H^^'(-E 2 + 2E 3 F jW ) 



(43) 



13 



Notice that 51 em / 5Ej = yields Ej = Fj xo and insertion of this equation back into I em 
recovers the original form of the Maxwell action. For the factor Fj 10 in the second term, 
we use the relation (I42h and express that in terms of the derivative with respect to r acting 
on Aj T . We then integrate this by parts to obtain 



Iem 



ttR [ drb-Wf^HWY-^ES 
J j 2X l 

7 3 



^3 

8vrG 5 

^3 

8ttG r 



(44) 



We note that the last term is a total derivative and does not contribute to the equations of 
motion. 

Finally, the scalar part of the action takes the form 



8vrG, 

^3 



8ttG 5 

Now in the total action 



2ttR j drbr 3 f- 1/2 H 2/3 Y^ 



1 AX 2 + R^Xi 



(45) 



I = I G + I GH + Iem + Is , (46) 

we clearly see that Aj T in Iem are the Lagrange multipliers. The variation with respect to 
those fields yields Gauss' law constraint of the system, 



b -l r S f l/2 H l/3j_ E 

J 2X 2 J 



0, 



(47) 



where the repeated index j is not summed over. We introduce the integration constants 



-iej so that 



E^^ibf-^H-^X 2 % 



(48) 



Iem 



(49) 



The constants ej describe the real (as opposed to complex) charges that satisfy Gauss' law 
and the factor of —% is inserted to account for the analytic continuation from Minkowski to 
Euclidean space. 

Using Equation ((351) . Iem can be expressed as 

U,R f dr6r- 3 r 1/2 ^- 1/3 E e i X l-^E(^)^} > 

3 3 

where the last term comes from the evaluation of the total derivative term in Equation 
and we have applied the boundary conditions (I28p and (|30l) . 
The total action at this point has the form 



^3 

8ttG r 



T=r B 



+ £k 2n I drbr3 r 1/2H2/3 { GT r-^G 5 r\) , 



(50) 
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where T T T is denned as 



epg l f{X' 3 Y _2_J_ 
r 6 # 4 M + R 2 X,-, 



(51) 



We have noted that G T T shown in (|33|) is free of the function b. Therefore, in the total 
action (|50p . we explicitly see that the function b is the Lagrange multiplier that enforces 
the constraint 

G T T = 8ttG 5 T t t . (52) 

One can verify that this, in fact, is the ( T r )-component of the Einstein equations and thus, 
this is the Hamiltonian constraint of the system ^ 



Given the expressions (|33p and (|5ip . the Hamiltonian constraint (|52p is a differential 
equation of functions /, H and Xj with respect to the variable r. We solve this differential 
equation by making ansatz which is similar to the one made in Reference |12j . 

2 2 

f(r) = 1 - J| + ^H(r) and Xj{r) = Hj 1 H{rf/\ (54) 

where tq is an undetermined parameter and the functions Hj are defined as Hj := 1 + 
qj/r 2 . Now because of the constraint (|15p imposed on the scalar fields, the function H(r) 
is determined to be 

H{r) = H 1 {r)H 2 (r)Hz{r) . (55) 

Thus we see that qj are the three parameters of the three fields; H and Xj with the constraint 
(|15p . Notice that this ansatz is consistent with the asymptotic condition (|22p . The form of 
the function H above implies that the topology condition ([19]) is equivalent to the equation 

f(r H ) = , (56) 

as advertised. Using this condition, we can trade the parameter ro of the function / with 
the parameter rjj, just as in Equation @. Then the Hamiltonian constraint (|52p is satisfied 
if and only if 

„2 



e] = qj{r H +qj] 



(57) 



9 In general, the Hamiltonian constraint of a gravitational system is given as G^n" — SirGsT^n 1 ' , 
where n u is a unit time-like vector. In our metric (|18[1 and action (|14fl , there is no non- vanishing off-diagonal 
components in the Einstein and energy-momentum tensors. Therefore, the Hamiltonian constraint of our 
system is G T T = 8nG^T T T which is the ( T T )-component of the Einstein equation. 

The variation of the action (|14|) with respect to the metric yields the Einstein equations of the system 

GW = Y l -^s[{d^Xj){d v X i ) - ^(dxXiX&Xi) - \g„uF jaP F^ + F ilt x F jvX \ + ■ (53) 

3 ' J 

From this equation, it is easy to verify that Equation (I52p with (|5ip is the ( T r )-component of the Einstein 
equations. 

The Einstein equations of the theory is frequently presented in the form that does not have the scalar 
curvature. Expressing the Einstein equations in this form requires to take the trace of the Einstein equation 
and this assumes the equality of other components than the constraint equation. Therefore, it is not adequate 
for extracting only the constraint equation. 
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where the repeated index j is not summed over. One observes that these relations are 
identical to the ones obtained in the on-shell analysis of Section [2] as shown in Equation . 

In solving the Hamiltonian constraint, we have adopted the rather tight ansatz (154j) . 
This narrow assumption on the functions restrict the form of possible solutions. Instead of 
adopting the ansatz discussed above, we can assume the forms of the scalars Xj and the 
function H only, and regard the Hamiltonian constraint as a first order differential equation 
of the function f(r). The solution to the differential equation generates a new family of 
solutions that includes the solution discussed in Section [2] as a special case. We will discuss 
this more general solution in Section 15.21 



3.4 The Boundary Reduced Action 



In the previous subsection, we have solved and eliminated the constraints of the system. 
Thus the action now takes the form 



8vrG 5 



f^H-^(r 3 H^)'] r=rB+ pY: {^)^ + H ! iy^ W ^} ■ (58) 



This action describes the thermodynamics of the black hole system in the box with the wall 
at r = r#. As a grand canonical ensemble, the temperature (5 and the chemical potentials 
jlj are defined at the wall and the grand canonical data at other locations in the box cannot 
be deduced from (3 and Uj, for the functions b and Aj T are the Lagrange multipliers and 
unspecified except on shell! 10 ! 

We are interested in the thermodynamics with the limit r# — > oo. The action has the 
asymptotic form 



SttGs I'rJ 



+ 5152 { E «? - 10 ('1« + <™ + <Ml)} + 0(4- 

i B 



U 3 



R 3 



8ttg s v*-< \r?) h 



(59) 



where we have defined 



E 



(60) 



with ro as appearing in Equation 

The second and third terms in the action are divergent as the wall approaches the 
boundary. Recall that we did not include the arbitrary background subtraction term in the 
Gibbons-Hawking term (I35p . This is because we prefer to regulate the divergence of the 



10 This can be physically understood. Imagine an experiment with a thermodynamic system in a box 
where the number of electrically charged particles of the system is allowed to vary. An experimenter sets 
up a temperature on the wall of the box using a kind of heat reservoir and also applies an electric voltage 
on the wall. The experimenter cannot define the temperature and the chemical potential of the system as a 
whole until the system reaches the thermodynamic equilibrium. The analogy is clear. 
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action by using the AdS / CFT- motivated counterterm technique [251 [251 [22] without arbi- 
trarily assuming a background geometry. Therefore, we employ the standard counterterms 



1 



8ttG 5 



d p 

dM -tC 



1 



dM 4 



8ttG 5 

where IZ4 is the scalar curvature with respect to the four-metric (f36l) and we have 

6 



r 2 B H(r B y/z ■ 



(61) 



(62) 



The evaluation of the counterterms is straightforward and yields the asymptotic form 

,4 

B_ 

'R? 



h+h 



^3 



re/ 



3^ + !( 1 +£§H + iE^ 



^ { E " 2 (<M2 + 9293 + g 3 gi)} + 



(63) 



We see that the first two terms in the above expression cancel the divergent terms in 
Equation (j59|) . 

In addition to the standard counterterms just computed, we will find that we need 
to include the "finite counterterm" discussed by Liu and Sabra [15], so that the resulting 
reduced action correctly reproduces the black hole solution of Section [2j Liu and Sabra 
argued that in the presence of matter fields, it is possible to consider other forms of the 
counterterms than the standard ones. The possible new counterterms include 



1 



nite 



where the scalar fields 



8vrG 5 j 9M 
1,02) are related to Xj [27], via 



and the vectors tij can be chosen as 



exp 



1 _ 

2 3 



a 1 



f-Lv^) , a 2 = (^=,-V2 



V^6' 



VV6' 



«3 



4 



(64) 



(65) 



(66) 



Under the assumption of the spherical symmetry, this is the only non- vanishing new coun- 
terterm. Though this term is finite as the limit, r B — ► 00, is taken, the inclusion of this term 
produces the mass consistent with the one reported in Reference [12]. To avoid diverting 
the topic here, we yield the further discussion on the finite counterterm to Section \5. 31 With 
the forms shown in Equations (|54p and (|55p . the finite counterterm (|64p is 



finite 



8ttG 5 \ P r B 



) ^2 { E <$ ~ (<M2 + + qsqi)} +0{^). (67) 



Combining all the results together, we have the action 
lu 3 R 3 \f 3 R\ 



I 



8vrG 5 



(68) 
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Recall from Section 13.21 that the inverse temperature (3 diverges as the limit rs — > oo is 
taken. The form of the action urges us to redefine the temperature as in Equation (|24p 
and similarly for the chemical potentials as in Equation (I29p . Adopting the definitions, the 
limit vb — ► oo yields the action 

I* =p{E-J2 - 2vr(^) 3 v^R) , (69) 

3 

where in the definition, the overall dimensionless constant W3-R 3 /8ttG^ has been absorbed 
[commonly in the literature, the convention G5 = (7r/4)i? 3 is adopted for the same effect] 
and the parameters qj and ej are assumed to be related through Equation (|57|) . We recall 
that the temperature j3 and jXj are not the physical quantities in the usual sense, and I* is 
an action that is defined only on the boundary. In this sense, this is the boundary reduced 
action derived from the bulk black hole thermodynamic system. 

We end this section by commenting on a few expected properties of the action. First, 
we note that extremizing this action with respect to rjj and qj, it correctly reproduces the 
chemical potentials ([8]) and the temperature (fl2|) of Section [2j Secondly, the constant 3/8R 
in the definition of E in Equation (|60p is the Casimir energy and it is known to exactly 
agree with the vacuum energy of the large N c N = 4 super- Yang-Mills theory defined on 
the manifold R x 5" 3 , through the AdS/CFT dictionary [25J. Apart from this important 
fact, this term does not contribute to the thermodynamics of the system. Other than this 
Casimir factor, the (off-shell) boundary reduced action has precisely the same form as the 
on-shell Gibbs free energy (|13j) with E being proportional to M in Equation (I10p Finally, 
the mean thermal quantities can be determined by taking derivatives of I* and evaluating 
them on shell yielding 

(E)=E, (Q j ) = e j /R 2 and S = 2n(^f yj H(r H ) ■ (70) 

The thermal quantities given above are functions of th and qj, and one can readily verify 
that the thermodynamic identity 

d(E) = ^dS + Y,Hd(Qj), (71) 
P 3 

is satisfied. 

4 The Global Structure 

We now examine the global behavior of the boundary reduced action I* in Equation (|69p . 
The action is a function over the parameter space of (rH,Qj) with the fixed parameters 
{(5,iij). The essential features are captured by the subspace where there is only one non- 
vanishing charge (and chemical potential), therefore, we concentrate our discussion here 
to this case. In Appendix |Al we show that the case with general charge configurations 
essentially reduce to the one charge case. In this section, we set 

R = l, (72) 

11 This remarkable result is also obtained in other systems such as the charged four dimensional black 
holes in asymptotically flat space [I], and in asymptotically AdS space [28]. Here we see that it is true even 
with the presence of the scalar fields that non-trivially couple to gravity and Maxwell fields. 
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and assume that all the length scales are measured in units of the (asymptotic) AdS curva- 
ture radius R. 



4.1 The Metastability 

We first present the main result; there is a metastable direction in the parameter space. 
Figure [2] is plotted for the configuration (ei,e2,es) = (e,0, 0) [or (91,42,93) = (9,0,0) 




Figure 2: The plot of the action with only one charge turned on. It is plotted against the variables ru and 
q, rather than e. The temperature and the (only one) chemical potential are set to (T, fj,) = (0.45, 1.05). The 
local minimum is located at (th , q) ~ (0.996, 1.23) and the unstable saddle point is at (rn , q) ~ (0.359, 5.38). 
The top flat part simply is a plot range cut-off. 

through the relations (J57J)]. The temperature T := 1//3 and the only chemical potential /i are 
set to (T, fi) = (0.45, 1.05). The action falls to minus infinity toward the direction (th, q) —> 
(0, 00), therefore, there is no true equilibrium state in this grand canonical ensemble and the 
local minimum at (r#, 9) ~ (0.996, 1.23) is the metastable black hole. There is an unstable 
saddle point at (th^q) ~ (0.359,5.38) and it is most likely that the metastable black hole 
decays by bouncing through this unstable state. In Figure [21 the action appears to go down 
linearly with respect to 9. In fact, the action has the rjj — ► limit 

/?{! + (! + 0(4). (73) 



I* 



This expression also shows that the global instability sets in at \i = 1. (Notice that Figure [2] 
is plotted for = 1.05 > 1.) 

The limit taken in Equation (|T3|) captures the essence of the metastability but it is 
slightly too simplified, so we now examine the metastable directions in details. For the one 
charge system in our current consideration, the action I* takes the form 



P 



{ 8 + 2 r ^ (1 + r H + l) +( l}- vArtir + ^ + ir) 



2ir r 



H 



H 



+ 9- (74) 
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We call the first term in the square bracket, "the mass part", the second term, "the EM 
energy" , and the last term in the above expression "the entropy" . As in the standard grand 
canonical ensemble, the mass part thermodynamically competes against the EM energy and 
the entropy. We are interested in the direction in the parameter space where the EM energy 
and the entropy overcome the mass part in the action. It is clear from the form of the action 
that this can possibly happen only in the direction q — ► oo. Expanding the action I* for 
large q, we obtain 

r=/3(l + \r\ - y/l + rjj p) q + 0(q l l 2 ) . (75) 

The next-to-leading order is q 1 ^ 2 but we concentrate on the leading term. For fixed values 
of /3, \x and q, this term has the minimum at m = as long as \x < 3, and for fi > 3, 
the minimum is located at rjj = %J /x 2 /9 — lF^I This means that the most likely direction 
of the instability is (rn,q) — > (0, oo) for 1 < n < 3 and is (rn,q) — > (\/^ 2 /9 — 1, oo) for 
fi > 3. Hence, it appears that the metastable black hole decays either by collapsing to 
rjj — > with diverging charge or into some finite horizon radius with diverging charge, 
depending on the value of /U. However, we suspect that this is not the physical behavior of 
the decaying black hole. First, the physical meaning of the value n = 3/R is highly unclear 
as oppose to /U = 1/R which has the relation to the BPS bound. The second reason is 
that our analysis is entirely classical and quantum corrections are not considered. Also we 
have ignored the effect of the back-reaction from the thermal gravitons. We should expect 
the modifications to the action due to those effects and should not trust the results too 
much in details without good reasons. Finally, our action can only describe the systems 
satisfying the assumptions that we have adopted along the derivation of the reduced action. 
Thus it is possible that a metastable black hole decays into some other spacetime with, 
say, a different topology or a geometry that is not spherically symmetric. We do know that 
there is a metastable direction in the parameter space but we cannot simply conclude the 
decaying process of the unstable black hole using the method adopted in this work. 

To summarize, our action shows that the metastability sets in at p = 1 as shown 
in Equation (I73h and a metastable black hole decays by accumulating large amount of 
charge due to the excessive value of the chemical potential. Though we have neglected the 
corrections to the action due to the quantum and back reaction effects, we expect that those 
effects do not modify the qualitative behavior of the system. In particular, as we discussed 
in the introduction, we have a good reason to expect the critical value of the chemical 
potential /i = 1 receives no corrections. It is commonly assumed that an instability of a 
system merely implies the existence of a new ground state but our analysis of the theory 
does not suggest such picture. 

4.2 Small Black Hole, Large Black Hole and Instabilities 

We now examine the behavior of the stable and unstable saddle points that the boundary 
reduced action possesses. For the one charge configuration, our action has a saddle point 
at the values of ru and q that satisfy 



g-!7 r ' + ', ,= ,M±^. (76) 
2r 2 H + 1 + q P V r% + q 



12 The coefficient of the next-to- leading order is —2-Krjj. Since the minimum is at th = as long as fi < 3, 
this order does not modify the critical metastability line fj, = 1 as in Equation (|73[) . 
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Since we are considering the grand canonical ensemble, these expressions should be solved for 
th and q in terms of given values of f3 and /x, and the solutions are the thermal expectation 
values of the equilibrium stateEl There are two relevant solutions with real and positive 
horizon radii. 

In Figure El rs of those solutions are plotted in the parameter space of (T, fi) . In the 




Figure 3: The two relevant solutions, rn(T, /i), of Equations (|76[) are plotted and superposed. The top 
leaf corresponds to the (locally) stable black hole (large black hole) and the other leaf to the unstable black 
hole (small black hole). The lower leaf has a cusp along the line of \x = 1 and the value of the horizon radius 
along this line is zero. The slight gap between the leaves is due to the numerical noise. 

blank regions of the parameter space, there are no saddle point solutions to the action. The 
blank region that surrounds the origin of the diagram corresponds to the part inside of the 
dotted line in the left diagram of Figure [TJ Also the blank region of very large fi is the local 
instability part of the phase diagram that corresponds to the region above the solid line of 
the left diagram in Figure [TJ Therefore, if we project Figure [3] down to the fi-T plane, the 
boundary of the plot is the (local) instability lines discuss in Section [2j The top "leaf" of 
Figure [3] is the locally stable black hole solution and the other bottom leaf is the unstable 
saddle point of the action. Those two saddle points of the action are commonly called "large 
black hole" and "small black hole", respectively, and we adopt this terminology here. In 
Figure HI we also plot the slices of the functions rji(T, //) and q(T,fi) against the chemical 
potential /i at the fixed temperature T = 0.45. 

Notice, in particular, the behavior of the small black hole at the line fj, = 1. At this 
critical metastability line, the small black hole, regardless of the temperature, shrinks to 
zero size in horizon and ceases to exist. Above the critical line, the small black hole gains 

13 The quantities should really be written as (rjj), etc., to be consistent with the notations adopted at the 
end of the previous section. However we simply use the notations th, etc., to avoid cluttering the equations. 
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Figure 4: Diagram (a) is plotted for the horizon radii th of the large and small black holes against fi 
at fixed temperature T — 0.45. The solid line corresponds to the large black hole solution and the dashed 
line to the small black hole. Similar plot is given in Diagram (b) for the charges q against fj, at the same 
temperature T — 0.45 as the other diagram. 



finite horizon radius again, but the charge jumps to a bigger value than that of the large 
black hole. The large and the small black hole saddle points of the action merge at the local 
instability lines (the solid and dotted lines in the left diagram of Figure [I]). 

In what follows, we describe the role of the small black hole in different phases. In 
Section [21 we explained that the region between the dotted and dashed lines in Figure [Q 
is the spinodal phase where the thermal AdS space without black hole is energetically 
preferred but the black hole solution remains locally stable. The action I* at the origin of 
the parameter space (th,q) corresponds to the thermal AdS space without black hole 1^1 
Therefore, in the spinodal phase, the large black hole saddle point does exist but the value 
of the action at this point is larger than the value at the origin. There is an unstable small 
black hole solution between the large black hole and the thermal AdS space in the (r^,g) 
parameter space. Thus, the large black hole in the spinodal phase is most likely to decay into 
the thermal AdS space, bouncing through the small black hole. The behavior is opposite 
in the high temperature side of the Hawking-Page phase transition line (the dashed line in 
Figure [1]) and below the metastability critical line fj, = 1. There, the black hole solution 
is energetically preferred to the thermal AdS space and the latter is most likely to decay 
into the large black hole, bouncing through the small black hole. In the region above the 
fj, = 1 line and below the local instability line (the solid line in Figure [1]) , the small black 
hole is located between the large black hole and the instability direction of the parameter 
space. Therefore in this region of the phase diagram, the large black hole is most likely 
to decay into the instability direction bouncing through the small black hole. The thermal 
AdS space in this phase can either directly decay into instability direction or decay through 
the large black hole. At the local instability line present above \i = 1, the large and small 
black holes merge and open up the instability direction without any barrier. 

To illustrate the described behavior of the saddle points, Figure [5] displays the sequence 
of the actions with respect to rising chemical potential and fixed temperature T = 0.5. 

14 Strictly speaking, the action does not describe the thermal AdS space because it has different topology 
from what is assumed in constructing the action. In fact, the origin of the action is not a saddle point, 
except when fj, = 1. However, it formally corresponds to the thermal AdS space with the correct equilibrium 
energy 3/8 and correct black hole parameters of th — = q. 
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The numerical data for the diagrams are summarized in Table [2j In the table, "LBH" and 
"SBH" stand for "Large Black Hole" and "Small Black Hole" , respectively. 



Diagram 


(a) 


(b) 


(c) 


(d) 




0.9 


1 


1.1 


1.4 


LBH (r H ,q) 


(1.208,0.7166) 


(1.211,1.000) 


(1.205,1.415) 


(-,-) 


SBH (r H ,q) 


(0.1087,0.1716) 


(0.000,-) 


(0.5048,6.885) 


(",-) 



Table 2: The numerical data for the plots in Figure [5] The abbreviations "LBH" and "SBH" stand for 
"Large Black Hole" and "Small Black Hole", respectively. 



(a) : n=0.9 (b) : fx=l 




Figure 5: The action is plotted for different values of the chemical potential. For all the diagrams, the 
temperature is set to T = 0.5 and other numerical data is shown in Table [2] The top flat parts are the plot 
range cut-offs. 

We comment that the decay processes discussed in this section are similar to the one 
analysed by Gross et al. [21] where their "Schwarzschild instanton" corresponds to the 
small black holes in our discussion^] The authors of Reference [21] discovered that the 
Schwarzschild instanton is an unstable saddle point by examining rather general pertur- 
bation to the black hole metric. Here we have restricted our consideration to the static 

The authors of Reference [21] also discuss the Jeans instability. But as they show in the reference, this 
instability is due to the "anti-screening" of gravitons that induces negative thermal mass squared. Therefore 
this mechanism is totally different from the one discussed in this section. 
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and spherically symmetric geometries. Moreover, we have made the very stringent form 
of metric ansatz and considered very small subspace of the full parameter space. As men- 
tioned near the end of Section [21 reducing the parameter space usually makes it hard to 
find the unstable directions. Therefore, it is rather remarkable to find the local and global 
instabilities in our narrow range of analysis. 

4.3 Below T 

In the regions of the phase diagram discussed above, temperature is always larger than Tq 
at which the horizon radius ru vanishes. The region of the phase diagram below To has 
been completely unknown because the analysis was done through the on-shell action which 
requires the existence of the black hole saddle point. Now that we have the off-shell action, 
we can discuss this region of the phase diagram. We have two results in this region. One is 
that the critical line at \x = 1 persists at all temperature including below To, and the other 
is that the black hole saddle points do not exist in this region. These two facts imply that 
below temperature Tq, a thermal AdS space decays directly into the instability direction 
without going through the black hole saddle points, provided that the chemical potential is 
larger than the inverse of the curvature radius. In fact, one can plot the action for a fixed 
temperature below To with chemical potential larger than one and finds a similar plot to 
Diagram (d) of Figure [5j Therefore, the critical line \i = 1 below Tq is not the metastability 
line of a black hole but it is the instability line of the thermal AdS space. 

5 Discussions 

5.1 AdS/CFT Boundary Temperature 

In Section 13.21 we found that the physical boundary temperature red-shifts to zero due 
to the asymptotic AdS geometry. Therefore, in defining the boundary action, we adopted 
rather ad hoc rescaled temperature given in Equation (|24"|) . We discuss this issue more in 
detail and argue that the rescaled temperature has a natural interpretation in the sense of 
AdS/CFT. 

The black hole temperature Tjj computed, for instance, by evaluating the surface gravity 
at the horizon or by Euclidean continuation to smooth geometry is not a physical temper- 
ature in general, and the physical temperature depends on the location of an observer. 
Therefore, the physical temperature (sometimes called Tolman temperature [29]) is given 
by 

T physica i = T H V(r B r 1 / 2 , (77) 

where we have assumed a spherically symmetric black hole geometry with radial coordinate 
r, red-shift factor V(r) and we locate the observer at r = rs- 

In the case of the asymptotically flat black hole spacetime, the red-shift factor V(r) 
can be given with the asymptotic behavior V —> 1 as r^^oo. Thus in this case, the 
Tolman temperature measured at the boundary of the space equals Tff and therefore this 
is commonly considered as the physical temperature on the boundary. The situation is 
different for the asymptotically AdS black holes. The red-shift factor in this case has the 
asymptotic behavior V — > r^/H 2 with R being the curvature radius of the asymptotic AdS 
geometry. This AdS red-shift factor causes the gravitational potential that prevents any 
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massive object from escaping to infinity and effectively acts as an infinite box [16J. Also we 
see in Equation (J77j) that the same AdS-term red-shift the physical temperature to zero at 
the boundary and thus Th is not the physical temperature of the boundary observer. 

We have mentioned in Section 13.21 that the AdS black hole thermodynamics in the in- 
finite volume has been discussed by adopting the "renormalized" temperature where the 
temperature is multiplied by the black hole red-shift factor so that the limit rs — > oo re- 
mains finite. However, since := lim^^oo Tphysicai^(^B) 1 ^ 2 = Tjj, the "renormalized" 
temperature is the temperature Th- The original work of Hawking and Page in Refer- 
ence [16| discusses the thermodynamic properties using (or Th) and also the following 
works by others have frequently adopted the same convention. (See, for example, Reference 
|17j and the references cited therein.) As emphasized by many, the redefined temperature 
Too is not the physical temperature of the boundary thermodynamics but as we stated in 
Section 13.21 there are indications that T^ somehow is the boundary temperature which is 
consistent with the thermodynamic properties. The reason why this is so has remained 
unclear. 

Here, we consider this issue in the light of AdS/CFT correspondence. Let us first recall 
how the conformal redundancy arises on the boundary of AdS space, based on Reference [8]. 
To make our discussion concrete, let us consider the AdSs space with the metric 

ds A dS = j^-yf [ - + y) 2dt2 + y'V + 4 ^°1] , (78) 

where the parameter R is a length scale of the geometry and the (dimensionless) time 
coordinate t is periodically identified. The boundary is given at y = 1 and the metric is 
ill-defined for this value of y. To obtain the metric that extends over to the boundary, one 
can multiply the AdS metric with a function r](x) 2 where the function rj(x) possesses a first 
order zero at y = 1. Again, for the sake of concreteness, we choose 

e u)(x) 

V{x) = —(l-y). (79) 

Here, oj is an arbitrary function over the spacetime including the boundary. Then the 
boundary metric ds 2 is given by restricting the extended metric to y = 1, i.e., 

ds 2 = [v{x) 2 ds 2 AdS ] = e 2uJ R 2 (-dt 2 + dQ 2 3 ) = e 2uJ (dr 2 + R 2 dn 2 3 ) , (80) 

where the function u now is understood to be defined over the boundary and the last 
equality is obtained through the analytic continuation to the Euclidean signature by setting 
r = iRt. Observe that the boundary metric is unique only up to conformal transformations 
due to the arbitrary function exp(w). We also see that the boundary topology is S 1 x S 3 . 
If we assume the existence of a conformal field theory on the boundary, the circumference 
of the S 1 time circle and the length scale of S 3 should be respectively interpreted as the 
inverse temperature and the length scale of the space on which the CFT is supposed to live. 
However, each of those quantities cannot have physical significance due to the conformal 
redundancy. A physical quantity of the field theory must be an invariant of the redundancy 
and in this case, it is given by the ratio of the S 1 circumference to the length scale of S 3 . 
In other words, the temperature of the boundary field theory should always be measured 
in the units of the length scale of the spatial three sphere. 
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The argument made in the previous paragraph closely parallels our problem of the black 
hole boundary temperature. The ill-defined nature of the AdS metric at the boundary is 
directly reflected to the fact that the local inverse temperature (3 defined in Equation (|23[) 
diverges at the boundary. Therefore, unlike the asymptotically flat case, the thermodynam- 
ics of AdS black holes is not well-defined in the same straightforward manner. In order to 
define the black hole temperature on the AdS boundary, one could extend the AdS metric 
as ds 2 in Equation (180j) and this is similar to the definition (124ft of the boundary black hole 
temperature (3. Notice that this could leave the definition of (3 ambiguous by an arbitrary 
function. However, the boundary metric has conformal symmetry and thus this black hole 
temperature cannot be a physical one. Therefore we realize that a possible physical black 
hole temperature at the boundary must be defined as the ratio between the circumference 
of S 1 time circle and the length scale of S 3 at the boundary; 

„ . ^-circumference at tb jH 

/^boundary := hm „ = lim = f3 . (81) 

r s ^oo length scale of b A at rs rg— >oo tb 

As we alluded in the last equality, the inverse black hole temperature [3 should really be 
interpreted in this way. [Equivalently, one may regard (3 as the rescaled temperature by an 
ambiguous function of the sort shown in Equation (|79|) . followed by taking the ratio of it 
to the length scale of S 3 , which is rescaled in the same way, at the boundary] 

Now the boundary temperature thus defined may not have any relation to the thermo- 
dynamics of the black hole in the bulk because it is not the Tolman temperature of the 
observer at the boundary!^ What one should be surprised here is that the boundary action 
I* with the temperature (3 does consistently describe the thermodynamics of the bulk AdS 
black hole. One could speculate from this that there would be a conformal theory on the 
boundary with the temperature (3 and this theory has identical thermodynamics as the bulk 
black hole, i.e., the AdS/CFT correspondence. 

The definition of the chemical potential in Equation (|29ft also is a conformally invariant 
quantity. To see this, notice that the integral of the gauge field one- form in Equation (|25ft is 
a conformal invariant. Then the product jlf3 is invariant under conformal transformations. 
Now since we have argued that the boundary temperature (3 is a conformal invariant, so 
must be the boundary chemical potential \x. 

5.2 More General Solution 

In Section 13.31 we solved the Hamiltonian constraint equation (|52p by adopting the ansatz 
shown in Equations (|54p and (|55p. Alternatively, we can take the ansatz for the functions 
Xj{r) and H(r) as before but do not assume the form of the function f(r). Then regard the 
constraint equation ([52]) as a first order differential equation of /(r) and get the solution 

r 2 H(r) 



3 r 4 + 2 r 2 £\ q . + qiq2 + q2q3 + qm 

3 3 



(82) 



16 We should comment that in Reference [30] . Witten has briefly mentioned the black hole temperature /3 
discussed here. However we are emphasizing that (3 is not a physical temperature of the black hole at the 
boundary in a usual sense. 
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where C is an integration constant and the length scale is set to R = 1. The constant C 
can be determined through the condition in Equation (|56|) and we have 



C = -E^4— -r 2 H U + 2j2<lj + 3r 2 H } . (83) 

3 3 

Unlike the previous solution to the Hamiltonian constraint equation, the parameters ej and 
qj are independent parameters of this solution. We can follow the steps taken in Section 
with the new function f|82j) and obtain the boundary reduced action 



I* =13 



J 3 H 3 



3=1 



2nr%^H(r H ) . (84) 



In extremizing this action with respect to the parameters (th, £j, Qj), one finds that there are 
one-parameter family of extrema. In fact, dl* /dm is directly proportional to ^ ■ dl* /dqj 
and the equations are not linearly independent. We can solve the six independent equations 
out of the seven dependent ones and express the solution, for example, as 

= 1 + X> + 2r " - ^H 3 {r H )- l yfW^), (85a) 

k^3 

e 3 = H( r H + Qj) ■ ( 85b ) 

Using a part of equations of motion, we can determine the Lagrange multipliers b{r) and 
Aj T (r) as 



ffWH- 1 '*, and AjT = -i^(^—-^—). (86) 

We thus have forms of all the fields and it is checked that they, together with the the relations 
(|85p . satisfy all the equations of motion including the Maxwell and scalar field equations 
with a free parameter, say, /30 Furthermore, we have computed the surface gravity at 
the horizon (with the same analytic continuation as the one mentioned in footnote 1 17|) and 
consistently obtained 

k = 2ir/p . (87) 

Also the thermodynamic identity can be shown to be satisfied, just as we have done at the 
end of Section 13.41 Parenthetically, the solution satisfies the smoothness condition (|2ip as 
we imposed in obtaining the reduced action. 

The asymptotic behavior of the function f(r) is 



/(r)=r 2 +(l + J>)+0(l 



17 We note that the gauge field shown above differs from Equation (0) by the factor of —if3/2ii which 
comes from the analytic continuation as shown in Equation (|27|) . 
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Up to this order, this agrees with the asymptotic behavior of the solution obtained in 
Section 13.31 and in particular, the first term is consistent with the asymptotic condition 
(|22p . One can check that the function f(r) takes the form in Equation (|54p if and only if 
the relations (|57j) hold. So while new f{r) is parametrized by the independent parameters 
(rn,ej,qj), the previous one is recovered by restricting the new function to the subset of 
the parameters by the relation ([57|) . As we have seen in Section 13.41 this relation reproduces 
the solution discussed in Section [2j Therefore, we have obtained the more general solution 
that include the previous solution as a special case. 

What is perhaps odd about this solution is that the grand canonical thermodynamic 
data, (f3,fij), cannot uniquely specify the properties of the black hole. Stated differently, 
the black hole can be uniquely characterized by the parameters, (E,ej,qi), where E is 
the thermodynamic mean energy derived from the action and we have chosen q\ to be the 
free parameter of the degenerate solution (|85p . Thus, the degeneracy seems to indicate 
a scalar hair of the black hole. But this conclusion sharply contradicts with the no-hair 
theorem of Sudarsky and Gonzalez [31J. They have shown that there cannot be scalar 
hairs for an asymptotically AdS black hole where the scalar potential approaches the global 
minimum near the boundary. Our scalar potential V in (|17p with the solutions Xj in 
(|54p actually has such asymptotic behavior. How can this be? One possibility is that the 
degenerate solution presented in this subsection is equivalent to the known solution through 
some reparametrization (we have not been able to find such reparametrization). Another 
possibility is that the no-hair theorem is somehow evaded in our theory, for example, due to 
the coupling of the scalar to the Abelian gauge fields (such coupling was not considered in 
Reference[31J). Clarifying this issue requires further considerations separate from the topic 
of this work. 

5.3 On Finite Counterterm 

In Section 13.41 we added the term called "finite counterterm" as shown in Equation (|64p . 
The original motivation for including the finite counterterm is as follows. Buchel and Pando 
Zayas |13j observed that when the mass of the -R-charged AdS black hole was computed 
using the standard counterterms, an unexpected nonlinear term with respect to the charge 
appeared. (The charge referred here is the parameter q and not e.) This nonlinear behavior 
of the charge in the mass can contradict the BPS-expected inequality between the mass and 
the charge, and more importantly, the thermodynamic identity is violated. By requiring 
the identity to hold, they added a term that consequently eliminated the nonlinear term. 
Then Liu and Sabra pointed out in Reference [H] that the inclusion of the counterterm 
associated with the presence of the scalar fields is possible. They showed that the mass 
computed with the finite counterterm (|64p reproduced the result obtained by Buchel and 
Pando Zayas. 

The counterterms that arise from the presence of matter fields are further studied by 
Batrachenko et al. in Reference [32] utilizing the Hamilton-Jacobi method. In this approach, 
the counterterms are treated uniformly throughout different spacetime dimensions and in 
particular, they have found that the scalar field counterterms are divergent, finite and 
vanishing for the spacetime dimensions four, five and higher, respectively. Note that the 
inclusion of the scalar counterterm in four dimension is necessary to render the action finite. 
Thus, it appears natural to include the scalar counterterm for the case of five dimensions 
but it still requires physical interpretation. 
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Liu and Sabra in Reference [15] mentioned that the option to include (or not to include) 
the finite counterterm in five dimension can be interpreted as the counter part to the choice 
of a renormalization scheme in a field theory. Since the inclusion of the finite countert- 
erm results in the mass that has the BPS-like mass-charge relation, they argued that the 
counterterm corresponds to the renormalization scheme that preserves the supersymmetry 
of the boundary field theory. If this interpretation is valid, the physical observables should 
not be affected by the finite counterterms just as the choice of the renormalization scheme 
does not affect the physical observables of a field theory. 

In our work, we have derived the action that describes the thermodynamics of the black 
hole observed at the boundary. The AdS/CFT conjecture postulates that the exponential of 
(minus) the action is the partition function of the dual field theory at strong coupling. Then 
the thermal mean values derived from the action, such as (E) or (Qj) of Equation (|70p must 
be the physical observables of the thermal field theory. However, we see that the addition 
of the finite counterterm does affect the thermal mean values. It is most clearly seen in the 
charges 




This expression is derived from the metric ansatz and the Hamiltonian constraint equation, 
so the form is independent of the counterterms. But the numerical value does change by 
the presence and absence of the finite counterterm because the location of the extremum is 
affected. Similar argument applies to the other mean values. Therefore, the observables of 
the theory do seem to depend on the finite counterterm. 

The dependence may imply that the finite counterterm has different physical interpre- 
tation. Another possibility is that the contributions to the partition function from other 
saddle points may make the physical observable independent of the finite counterterm. In 
any case, the observation made here urges us to investigate the physical interpretation of 
the finite counterterm and better understanding is desired. 

5.4 Comparison to the Field Theory Analysis 

In this work, we have found an additional structure in the phase diagram of i?-charged black 
hole system. The phase diagram in Figured] is now augmented by the line \i = 1/R that 
separates the stable and metastable black hole states for T > Tq and below To, the same line 
is the instability critical line of the thermal AdS space. As mentioned in the introduction, the 
existence of the metastable black hole has been anticipated from the AdS/CFT-conjectured 
field theory analysis in Reference [9] . The current status of the phase diagrams in the gravity 
and field theory sides is schematically summarized in Figure [6] We see the agreement in 
the general structure of those phase diagrams and the additional match in the metastability 
critical line fi = 1/R is the improvement made in this work. 

In Section 14.31 we saw that in the region of the phase diagram below To, a thermal 
AdS space decays directly into unstable direction without going through the black hole 
saddle points. The corresponding field theory analysis has not been carried out at the low 
temperature. However, given the result of the gravity analysis, we have some expectations 
for the behavior of the field theory. First we note that the "confinement/deconfinement" 
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Figure 6: The schematic phase diagram of the _R-charged black hole is shown on the left-hand side and on 
the right-hand side, that of the field theory which is supposed to be the AdS/CFT dual is shown. According 
to the AdS/CFT correspondence, the black hole system is identical to the strongly coupled region of the 
field theory, as indicated in the diagram. For both figures, jj, is the largest of the three chemical potentials. 
In the strong coupling diagram, To is the temperature at which the black hole radius reaches to zero. Its 
value depends on the pattern of chemical potentials and equals (irR)~ 1 for a single non-zero charge, while 
To — for three equal charges. On the other hand, the nature of the instability in the low temperature 
regime of the weakly coupled field theory is yet to be explored and is indicated by the dashed line in the 
diagram. 



phase transition line is obtained in the zero coupling limit of the field theory! 8 I The phase 
transition line and the \x = 1 line merge exactly at T = 0. The zero coupling limit is 
taken so that the resulting free theory is continuously connected to the theory with finite 
coupling constant 1^1 Thus we expect the general structure of the phase transition line to 
remain the same at a finite coupling, but it is likely that the details would change. For 
instance, it is possible that the phase transition line at the finite coupling is modified so 
that it merges with the instability line at some finite temperature. This behavior could be 
similar to the one observed in the gravity analysis. To verify this expectation, one needs to 
find the confinement /deconfinement phase transition line at weak but finite coupling with 
nonzero chemical potentials and also the instability line at the low temperature region must 
be determined. We believe that the analysis, especially in the low temperature limit, is 
doable and the work is in progress. 

In the field theory analysis of Reference [9], it has been found that the decay of a 
metastable state is most likely to occur by one of the eigenvalues of the scalar field splitting 
up from the rest of the eigenvalues in the field space and indefinitely growing. Our analysis 
in this paper does not necessarily suggest this mechanism. If we attempt to interpret this 
phenomenon in terms of string theory, it would correspond to a D3-brane separating from 
a stack of rotating D3-branes due to the high angular velocity of the D3-branes in their 
transverse directions. Therefore, we could also detect the metastability by examining the 
behavior of a probe D3-brane immersed in the background of the i?-charged black hole or 
of the rotating D3-branes (before taking the decouping limit). This scenario is currently 

18 The N = 4 super- Yang-Mills theory does not confine. The terminology "confinement" is used in accor- 
dance with the behavior of an order parameter, the Polyakov loop. See Reference [9] for more explanations. 

19 Technically, the zero coupling limit is taken while keeping the Gauss law constraint. 
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being investigated. 

Both gravitational and field theory analysis show that two theories cease to exist above 
the critical value of the chemical potential where the local instability sets in. In the 
AdS/CFT context, this means that the theory at strong and weak coupling regimes both 
become unstable at a large value of the chemical potential. One might naively interpolate 
the results and speculate that the theory in any coupling regime is unstable beyond a cer- 
tain value of the chemical potential. Though we think that this is the most likely scenario, 
we still cannot completely exclude the possibility that the theory develops a new ground 
state beyond the instability line and defines a Higgs phase. As mentioned in Section [H 
we did not take into account of the back reaction from the thermal radiation and also the 
quantum corrections are neglected. (Needless to say, stringy effects have not been consid- 
ered which we expect to be important when th becomes comparable to the length scale of 
string theory.) On the field theory side, the analysis has been carried out at one- loop level 
and higher loop with finite N c effects are still unknown. They may or may not change the 
global behavior of the potential in the theory. 
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A Metastability in General Charge Configuration 

In the main text, the metastability of the i?-charged black hole is examined in the subspace 
of the full parameter space with the only one nonzero charge. This appendix augments the 
argument that the metastable black hole system with a general charge configuration tends 
to have the unstable direction into the single charge configuration. 

As in the single charge case of the main text, a possible instability occurs only when 
one of the charges diverges. This is because the tq factor in E of Equation (j69H overwhelms 
the rest of the terms when more than one charges become large. Therefore, without loss of 
generality, we expand the action with respect to large q\ to obtain 

r=fi(l+ Z -x\ - ^/l + xl Ml ) qi + 0(q\ 12 ) , (90) 

where we have defined the parameter x\ := rf^Hi^r^H^^n)- This is exactly the same 
form as in Equation (|75p and therefore, for fixed values of (3, [i\ and qi, the leading term 
has the minimum at x\ = as long as /ii < 3, and for [i\ > 3, the minimum is located 
at x\ = -\J Hi/9 — 1. This implies that as long as 1 < [i\ < 3, the unstable direction is 
( r H, Qi, <Z2) Q3) — >■ (0, 00, 0, 0) and for [i\ > 3, the parameters rj{, qi and g3 may be nonzero 
but none of them may diverge. One of the consequences of this consideration is that the 
constrained system with qi = qi , q^ = or q\ = qi = 53 does not have the global instability 
(and as mentioned near the end of Section [21 these configurations do not have the local 
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instabilities as well). Those constraints reduce the dimension of the parameter space and 
the unstable directions are closed. The discussion made in this paragraph is illustrated in 
Figure [7] for the case with the two nonzero charges. Observe that the instability direction 

1.5 



1 

0.75 

r 0.5 

0.25 



q\ 2 

Figure 7: In this plot of the action, the configuration with qz = is chosen. The temperature, the horizon 
radius and the chemical potentials are set to T = 0.6, rn = 0.01 and (/ii,/H2) = (1.2,0.85), respectively. 
The diagram shows that an instability exists only in the direction of the large charge whose corresponding 
chemical potential exceeds the critical value 1. It also exhibits that the direction along q\ — q2 is stable. 
The top flat part is a plot range cut-off. 

opens up for the large charge whose conjugate chemical potential is larger than one. In 
this plot, it is the direction of large q\ and not q<i- Also note that the action does not have 
instability along the line q\ = qi- 

The result obtained in this appendix confirms the claim made in the main text; the 
metastability in the case with general charge configurations essentially boils down to the 
system with one charge. 
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